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Abstract 



A method is presented for the computation of the one-loop effective action 
at finite temperature and density. The method is based on an expansion in 
the number of spatial covariant derivatives. It applies to general background 
field configurations with arbitrary internal symmetry group and space-time 
dependence. Full invariance under small and large gauge transformations is 
preserved without assuming stationary or Abelian fields nor fixing the gauge. 
The method is applied to the computation of the effective action of spin zero 
particles in 2+1 dimensions at finite temperature and density and in presence 
of background gauge fields. The calculation is carried out through second 
order in the number of spatial covariant derivatives. Some limiting cases are 
worked out. 
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I. INTRODUCTION 



In their pioneering work, Deser, Jackiw and Templeton Pj] noted that gauge theories in 
odd-dimensional spaces naturally admit a local term of topological nature, known as Chern- 
Simons term || (see || for a recent review). One of the interesting properties of the non- 
Abelian Chern-Simons term is that under gauge transformations it changes proportionally to 
the winding number of the transformation. Thus, when the action contains a Chern-Simons 
term, the partition functional of the system is only well defined if the coupling constant of 
the Chern-Simons term is properly quantized. 

It was later realized [[§-[§] that the Chern-Simons term is induced by quantum fluctuations 
when gauge fields are coupled to odd- dimensional fermions. Such term comes out with the 
correctly quantized coupling constant and so full gauge invariance is preserved (although 
possibly at the price of spoiling parity invariance Because the Chern-Simons term is a 
polynomial in the gauge fields and their derivatives, it can be obtained through a combination 
of perturbative and derivative expansions. 

In the so-called imaginary time formalism for field theory at finite temperature |7|,|8[], the 
space-time has a non-trivial topology, since the time is effectively compactified to a circle. 
This allows the existence of topologically large gauge transformations even in the Abelian 
case. When the problem of the induced Chern-Simons term is studied using the method 
just mentioned of retaining a low number of fields and of derivatives, a puzzling situation 
appears, namely, the coefficient of the Chern-Simons term turns out to be a smooth function 
of the temperature, and hence it violates the quantization condition ||. The situation has 
recently been clarified by considering a simple + 1-dimensional model [|Hj which can be 
computed in closed form. There it is seen that full gauge invariance holds for the exact result 
but it is broken by perturbation theory. This is not difficult to understand, since in simple 
cases gauge invariance under large gauge transformations is equivalent to periodicity of the 
effective action as a function of the gauge field, whereas perturbation theory corresponds to a 
Taylor expansion of that function. Clearly, the property of being periodic is not maintained 
in general by a truncated Taylor expansion. In [] 1 it was noted that full gauge invariance 
is always a property of the exact result, since it follows straightforwardly from using a (- 
function regularization. In ]T2|] the exact effective action of fermions in 2 + 1 dimensions was 
obtained for the case of Abelian and stationary background gauge fields. 

The problem of preserving full gauge invariance at finite temperature is not tied to 



odd-dimensional theories [13-17| nor to fermions [18,19]. It appears whenever perturbation 



theory is involved. This is unfortunate, since, as noted in ||, "at finite temperature, per- 
turbation theory is one of the few tools we have." In this work we show that it is possible 
to carry out detailed calculations of the effective action fully preserving gauge invariance, 
without restricting oneself to particular configurations such as Abelian or stationary ones, 
and without choosing a particular gauge. The study of simple cases [|iT)|-|T2]] shows that 
the problem with gauge invariance comes through the scalar potential A (x). The finite 
temperature effective action is non-local in time but it is local in the space variables. This 
suggests to consider an expansion in the number of spatial covariant derivatives only. The 
time component is treated non-perturbatively in order to avoid destroying gauge invariance. 
(See [f2(J for another discussion of derivative expansions at finite temperature.) 
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It should be emphasized that the expansion in the number of spatial covariant derivatives 
is not tied to a particular method of computation, since it can be obtained from the exact 
result (namely, by considering an appropriate spatial dilatation of the background fields) 
and is fully gauge invariant. What it is shown here is that it is also amenable to explicit 
computation order by order, through a combination of the method of symbols and ^-function, 
for instance. This combination works very well and has been applied at zero temperature 
for fermions with local GTI] and non-local actions At finite temperature it has been 



applied to fermions in odd-dimensions |23j as well as in even-dimensions ||17|| , however, for 
technical reasons, this has been done choosing a particular gauge. In the present work we 
remove the necessity of any choice of gauge. It turns out that previous formulas in |23| . |T7| 
can be reinterpreted and rewritten in a manifest gauge invariant form. 

Although a naive perturbative expansion in Aq breaks gauge invariance, the very method 
of calculation suggests an expansion in powers of the temporal covariant derivative in the 
adjoint representation which preserves gauge invariance. This yields the remarkable result 
that, even at finite temperature, the theory is local if expressed in the appropriate variables. 

The method is explicitly applied to the case of relativistic scalar particles in 2 + 1 di- 
mensions. The computation is carried out through second order in the number of spatial 
covariant derivatives. Contact is made with the relativistic Bose gas. 



II. GAUGE INVARIANT DERIVATIVE EXPANSION AT FINITE 

TEMPERATURE 

A. The mathematical problem 

The aim of this section is to present a scheme to address the computation of the one-loop 
effective action at finite temperature preserving gauge invariance at every step. (We note 
that gauge invariance in this work refers always to vector gauge transformations.) To fix 
ideas consider the case of scalar particles in d + 1 dimensions in presence of background 
gauge fields. This case will be worked out later for d = 2. The Euclidean action of the 
system is 

s = J d d+1 x ((zv>) f (/v) + ™V</>) (i) 

where = <9 M + is the covariant derivative. The finite temperature condition can 
be implemented by using the imaginary time formalism, that is, by compactification of the 
Euclidean time to a circle so that the fields <fi and A^ are periodic functions of xq with period 
(3 = 1/T (T being the temperature). After functional integration over <p(z), the Euclidean 
effective action is formally given by 

W s [m, A] = Tr 6 log(-£>J + m 2 ) . (2) 

The subindex b recalls that the functional trace is to be taken in the Hilbert space of bosonic 
wave functions, i.e. with periodic boundary conditions. 

Presently the mathematical problem to be addressed is the computation of quantities of 
the form 
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F[M,A]=Tr(f(M,D)) 



(3) 



where is the covariant derivative and M(x) collectively denote one or more matrix 
valued functions of x M representing other external fields in addition to the gauge fields. The 
trace refers to the Hilbert space 7i of wavefunctions with space-time and internal degrees of 
freedom, the space-time manifold has topology AAd+i = S 1 x Ai d and the wavefunctions are 
periodic for bosons and antiperiodic for fermions. 

Although slightly pedantic, it will occasionally be convenient to regard M not as functions 
but rather as multiplicative operators in Ti, i.e., operators commuting with the operators 

and otherwise with arbitrary structure in internal space. Likewise are differential 
operators of the form + with A^ multiplicative. The quantity f(M, D) denotes an 
operator constructed out of M and in the algebraic sense, that is, f(M,D) is a linear 
combination (or series) of products of M and multiplied in any order with constant 
c-number coefficients. In order for M and to be well-defined operators in Ti, M(x) and 
A^{x) are required to be periodic functions of x . In addition we will assume that the fields 
are sufficiently convergent at infinity and the function / is well-behaved. This means, in 
particular, that / is one-valued and sufficiently convergent at infinity as a function of to 
ensure the existence of the trace (by avoiding ultraviolet divergences). 

A gauge transformed configuration (M u , A u ) is one of the form 

M u (x) = U~ 1 (x)M(x)U(x) , A"(x) = U~ 1 (x)d fl U(x) + U- 1 (x)A IM (x)U(x) , (4) 

where the gauge transformation U(x) is a periodic function of x$ which takes values on 
matrices in internal space. This corresponds to a similarity transformation of D^, namely, 
= dp + A^ = U^DyU where U is to be regarded as a multiplicative operator in 7i. 
Because f(M,D) is constructed with M, D and c-numbers, it follows that f(M,D) also 
transforms under a similarity transformation 

f(M u ,D u ) = U- 1 f(M,D)U (5) 

and so 

T[M U ,A U } = T[M,A] (6) 

using the cyclic property of the trace, which holds due to our regularity assumptions for 
f(M,D)fl 



1 In practice, T[M, A] is only computed for a subset of configurations (M, A) and only the subgroup 
of gauge transformations which leave invariant such a subset are relevant. For fermions, the internal 
space includes Dirac space as well flavor degrees of freedom, and the 7 M matrices are included in 
M (they are not c-numbers). In this case only gauge transformations in flavor space are relevant 
since they are the ones that preserve the form of 7„. 
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B. The method of symbols 



Assuming that the operator / = f(M, D) admits a complete set of eigenfunctions, f\n) = 
A n |n), the functional trace is simply T[M, A] = ^2 n X n - In this form gauge invariance is 
obvious since / and f u are related by a similarity transformations and hence they have the 
same spectrum. 

The gauge invariance of T[M, A] is also manifest computing the trace in the basis \x) 
of eigenfunctions of x^, normalized ) = 8(x — x') (a periodic delta function in the 

temporal direction) 

T[M,A] = J d d+1 xtr(x\f(M,D)\x) (7) 

(where tr refers to internal space) because U is a multiplicative operator and so tr(x\f(M, D) \x) 
is gauge invariant without integration over x. However, computationally it is more conve- 
nient to use a basis in momentum space \p) 

(x\p) = e p * , (p\p') = f36 pop ,(27r) d 5(p - p') , (8) 

(to avoid unessential factors of i we take the convention of using purely imaginary momenta 
p M but J d d p below denotes the usual integral in R d and 6(p — p') denotes the corresponding 
delta function). The frequency takes the Matsubara values po = 2nin/f3 for bosons and 
Po = 2ni(n + \)/ (3 for fermions. Note that we have assumed that the space manifold Add 
has a topology IR d . In this basis 

r[M, A] = ^E/ -^*W(M,D))p). (9) 

At this point the symbols method can be used (see e.g. f21| , |2"4|| ): let |0) denote the 
state with p = 0, then using the identities \p) = e xp \0) (where e xp acts as a multiplicative 
operator and the quantities p M are constant c- numbers) as well as e~ a;p Z) At e xp = _D M + p^, 
e~ xp Me xp = M, one obtains 

(p\f(M,D)\p) = (0\e- xp f(M,D)e* p \0) = (0\f(M, D + p)\0) , (10) 

and so the functional trace can be cast in the form 

r [M,A] = l]T I ^- d U(0\f(M,D + p)\0). (11) 

In this expression it is clear the requirement of regularity on /: the functional trace comes 
after integration over momenta and sum over frequencies and this requires / to be sufficiently 
convergent for large p M . Let us remark that |0) is periodic rather than antiperiodic in the 
temporal direction. The information on whether we are dealing with bosons or fermions is 
now contained solely in the values taken by po- The state |0) satisfies 

(x|0) = l, ^|0) = (0|^ = 0, (0|0) = Jd d+l x. (12) 
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In addition, when h is a multiplicative operator, h\x) = h(x)\x), 

(0\h\0) = [ d d+1 xh{x). (13) 



It follows that appearing inside f(M,D + p) in eq. ( |TTD acts derivating everything to 
its right (or its left, by parts) and then vanishes after it reaches |0) (or (0|). This is a 
well-defined working rule and from it stems the usefulness of the symbols method. 

Unfortunately, gauge invariance is no longer manifest when using the momentum basis. 
In fact, tr(0|/(M, D + p)|0) is not gauge invariant because |0) (or more generally \p)) is not 
covariant under local transformations. For instance, according to the rules given in eq. (O) , 



tr<0|Dj|0> = tr(0|Aj|0> = J d d+1 xtr[Al(x)] (14) 

breaks gauge invariance. However, 

tr(0|[Z> M> A,] 2 |0) = J d d+l xtr[F^(x)] (15) 

does not. Note that [D^, D u ] is a multiplicative operator whereas D 2 is not. As a rule, 
when an operator g(M, D) (a gauge covariant operator) is multiplicative, eq. flT3"|) applies 



and tr(0|g(M, D)\0) is gauge invariant |^TJ. In eq. ( [TI]) gauge invariance is only recovered 
after integration over momenta and sum over frequencies.^ This will be further discussed 
subsequently. 



C. The derivative expansion at finite temperature 

By computing the functional trace we essentially mean to end up with purely multiplica- 
tive operators, since this implies that the functional is expressed as the integral of a function 
over space-time. At zero temperature this is usually equivalent to saying that all derivative 
operators appear inside commutators. In addition it means to carry out as many implied 
sums and integrations (over frequencies and momenta or other parameters) as possible. 

In general it is not possible to compute T[M, A] in closed form and one must resort to ap- 
proximations. The standard approach is to make power expansions in one or more operators 
appearing in f(M, D) while the remaining operators are treated non perturbatively. As will 
be clear below, a naive expansion in powers of D would break gauge invariance, therefore, 
because it is in general difficult to work with two or more non perturbative operators un- 
less they are commuting, and our present emphasis is in the preservation of manifest gauge 
invariance rather than in a particular computation, we will keep Dq as the only operator to 
be treated non perturbatively, and expand in all other operators M and D. 



2 An elegant method has been presented in [M which yields gauge invariant expressions prior to 
momentum integration, at the price of introducing derivatives with respect to The method has 
not yet been extended to include discrete momenta, as required at finite temperature but in can 
be applied to the integration over p. 
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Before proceeding, let us be more precise about the meaning of expanding in powers of 
M and D. A convenient way to define the expansion is by introducing constant c-number 
bookkeeping parameters, f(M, D , D) — > /(A]M, D , \ 2 D), so that counting powers of those 
operators is equivalent to counting powers of Ai i2 . This procedure preserves gauge invariance 
since it amounts to a modification of the function /. After applying the symbols method 
(through — ► D^ + p^, cf. eq. ([□])) the factor X 2 will affect D and p, however, this can be 
brought to the form /(AiM, Dq+pq, X 2 D+p) by a redefinition of p. Therefore the expansion 
can be formulated as an expansion in powers of M and D in f(M,D +p)-f\ Equivalently, 
the expansion in D can be obtained directly from the functional T[M, A] by means of a 
covariant spatial dilatation, namely, M(xq,x) — > M(x , X 2 x), A (x ,x) — ► A (x , X 2 x) and 
A(xq, x) — > \ 2 A(xo, \ 2 x). This guarantees that the expansion is well-defined, i.e., it depends 
on the functional itself and not on how it is written or computed. 

The situation is completely different for an expansion in powers of D . A bookkeeping 
parameter D — > X 3 D can be introduced in f(M,D) and this defines a new A 3 -dependent 
gauge invariant functional. Nevertheless this functional is not useful since it presents an 
essential singularity at A3 = 0, as can be seen in the simple case of fermions in 0+1 di- 
mensions |1(J. (The dependence on A1.2 is analytic or at least asymptotic under suitable 
regularity conditions on the fields and on the function /.) After applying the symbols 
method, f(M, X 3 D + ^3Po ; D) is obtained. However, because p is a discrete variable this is 
not equivalent to f(M, X 3 D +p , D). Therefore expanding in the explicit D in f(M, D+p) 
does not correspond to a modification of / and in fact violates gauge invariance. Also, it is 
not possible to introduce A3 by means of a rescaling of type xo — > A3X0 of the field configura- 
tion (M, A) since this transformation violates the periodicity condition on the wavefunctions 

of n. 

After expansion there will be all kind of terms which will be products of single factors of 
M and D as well as operators depending non perturbatively on Dq (by this we merely mean 
that all orders of Dq are retained). It is always possible to bring all D operators to the right 
producing commutators, so that we end up with two kind of terms: i) terms in which all 
operators D appear only in commutators (more precisely, in the form [D, ]) and ii) terms 
with unsaturated factors D at the right (i.e., D not inside a commutator). The terms of the 
first type are multiplicative operators regarding as-space, although they are still differential 
(or pseudo-differential) operators with respect to xo-space. The terms of the second type 
are non-multiplicative in a>space. As we have argued above, these latter terms break gauge 
invariance and in fact they will cancel after integration over p. This can be seen as follows: 
let us replace D by D + a, where a is a constant c-number. This replacement has no effect 
on the terms where D is in commutators, but counts the contribution from the terms with 
unsaturated D. However, it is clear that there is no such a contribution after integration 
over momenta since a can be compensated by a similar shift in the integration variable p. 
Thus at the end, all operators D appear in commutators only. (The same result is obtained 



3 Of course, the actual expansion in powers of D must be done after arriving at eq. ( |il| ) (or other 
similar formulas in other approaches, such as Schwinger proper time method) i.e., after has 
become + p^, since otherwise powers of p^ would be generated as well and that would destroy 
the ultraviolet convergence of the formula. 
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directly using the method of Pletnev and Banin p4[ .) A similar argument would break down 



for D : a shift to -Do + «o cannot, in general, be compensated by a shift in p since at finite 
temperature the frequency is a discrete variable. 

From the previous discussion it follows that we only have to retain those terms where all 
operators D are in commutators. Because all operators are now multiplicative in a?-space, 
x becomes just a parameter in what follows. Let us consider a typical term: 

TT = i^y -^tx{0\a 1 (D +po,p)Xa 2 (D +po,p)Ya 3 (D + p 0l p)\0). (16) 

X and Y are multiplicative and gauge covariant operators constructed with and M, 
and the ai(x,y) are some functions. At this point the integration over p M is non trivial 
(even for the simplest forms of the functions aa(x,y)) because p M appears in different and 
non-commuting operators. A possible approach is to express the operators in terms of their 
matrix elements using as basis a complete set of eigenstates of D (in the Hilbert space of 
time and internal degrees of freedom). These matrix elements are then ordinary functions 
of pp. Instead of that, we will use the equivalent prescription of labelling the operators Dq 
according to their position with respect to X and Y: the symbols Dqi, Dq 2 and Dq 3 will be 
used to denote the operator D Q is positions 1 (before X), 2 (between X and Y) and 3 (after 
Y) respectively. In this notation 

TT = i^ J -^tT(0\a 1 (D 01 +p ,p)a 2 (D Q2 +p Q ,p)a 3 (D 03 +p ,p)XY\0). (17) 

An immediate consequence is that the labeled operators are effectively commuting and 
the momentum integration and frequency summation can be carried out as for ordinary 
functions. The result can be written as 

TT = tr((%(A)i, A)2, D 03 )XY |0> (18) 

where the function g is defined by 

1 x f d d p 

9{x,y,z) = J j^yai{x+Po,p)a 2 {y + po,p)a 3 {z + po,p). (19) 

(Note that there will be two versions of g, the bosonic one and the fermionic one, which are 
related by a shift of m / ' f3 in their arguments.) By construction the function g is periodic: 

2ni 2ni 2iri 
g{x,y,z) = g(x + —,y + —,z+—). (20) 

This is an immediate consequence of the sum over Matsubara frequencies and ultraviolet 
convergence of the expressions. 

This periodicity property is essential to codify the gauge invariance of the original ex- 
pression. To see this, let us introduce the operation T>^ which is defined as T>^X = [D^jX] 
for any operator X. Consistently with our previous notation, we will denote by Vqi the 
action of T>q in position 1 (i.e., on X) and by T>q 2 the action of T>q in position 2 (i.e., on Y). 
Then clearly 



S 



(21) 



The interesting point is that these formulas hold for arbitrary functions of V m and V 02 as 
well. This follows from the well-known identity e A Be~ A = e' A ' ^B: for any c- number A 

e x(D (n -D 02 ) XY = e \D OXe -\Do Y = (e^X) Y = e XT>01 XY , (22) 

and this identity immediately extends to arbitrary functions of _D i — D 02 , and analogously 
for V 02 . This allows to make everywhere the replacements 

D 02 = An - V 01 , D 03 = An - An - V 02 (23) 

and use D 01 , V 01 and V 02 as the independent variables to work with. The advantage of 
doing this is that the action of "Doi an d T> 02 on X and Y produces multiplicative and gauge 
covariant operators. On the other hand, the presence of the operator Dqi, which is outside 
commutators, combined with the gauge non-covariant operation (0| |0), still can introduce 
gauge non-invariant contributions. This is avoided thanks to the periodicity property of 
g(x,y,z) as we will show now. Indeed, the periodicity property allows to write the term as 

TT = tr(O\ip{e- 0Do \V ol ,V O2 )XY\O) , (24) 

where the function ip is defined by 

<p(e~ l3x ,y,z)=g(x,x-y,x-y-z). (25) 

The periodicity condition of g ensures that the function tp(uj,y,z) is one-valued (it depends 
on uj and not just on log (a;)). In order to bring the expression into a manifestly gauge 
invariant form, we will use the following property: 

e -tD = e -W Te - Q +t A ( X ' ,XW ^ t > Q (26) 

Here t is just a parameter and T denotes time ordered product. (The quantities D , do and 
represent operators and the product refers to a product of operators so that do is not directly 
derivating x .) This equation can be easily proved by the standard procedure of showing 
that the two expressions satisfy the same first order differential equation in t and coincide 
at t — 0. The left-hand side is a manifestly gauge covariant operator. It is interesting to see 
how gauge covariance is realized in the right-hand side: the time-ordered product from xq 
to x + 1 (with x fixed) transforms with U (x , x) at the right and U^ 1 (x + 1, x) at the left, 
and this latter factor is transformed into U^ 1 (xq, x) after commutation with e~ t9 °, thus the 
product of the two factors transforms covariantly at (xq,x), as e~ tL>0 . 

In particular, by taking t — (5 in the previous formula, one obtains the identity 



where 

"Xo+13 



/ rxo+P \ 
Q(x) = Texp I - / A (x' , x)dx' J . (28) 



9 



(Again e _/39 ° Q is to be understood as the product of two operators.) Beyond the interval 
[0,(3] A (x) is defined as a periodic function of the time, so Q(x) is also periodic. Although 
Q(x) is non-local in terms of A , it behaves as a local field which takes values on the gauge 
group. In particular it transforms covariantly at x: 

tf(x) = U-\x)n{x)U{x) . (29) 

The matrices Q(x) at different values of Xq, but equal x, are related by similarity trans- 
formations and their trace, the Polyakov loop, is independent of Xq. Another important 
property is 

V Vl = [D Q , Q] = . (30) 

On the other hand, the effect of the operator exp(— (3 do) is to produce the shift xq — > 
Xo — (3, therefore it is equivalent to the identity operator on the space of periodic functions in 
which we are working (as noted the periodic wavefunction |0) appears regardless of whether 
we are considering bosons of fermions). So in this space 

e^ D ° = Q . (31) 

This produces the manifestly gauge invariant expression 

TT = tr^^Doi, A) 2 )Ar|0> . (32) 

(The label 1 in Q indicates to put this operator in position 1. The relative order between 
Voi and fli is immaterial due to eq. (|30|) .) 

It can be noted that all previous manipulations, starting from eq. ([TED , hold also without 
taking tr(0| |0). Inside tr(0| |0) integration by parts implies that Z>oi is equivalent to — V 02 
(or equivalently, that D 03 = D i), and so we have the final formula 

TT = ti(Q\<p{n 1 ,V O2 )XY\0) , (33) 

where 



( p{<u,y) = (p{w,-y,y). (34) 

The whole point of these manipulations was to end up with a manifestly gauge covariant 
and multiplicative operator so that tr[<^(fii, Vo 2 )XY] is just a gauge invariant function of 
x, constructed with the fields A^{x) and M(x): 

TT = J d d+1 xtr[(p(tt 1 ,Vo 2 )XY] , (35) 

The fact that g(x,y,z) is periodic is essential to produce the gauge covariant operator 
Q(x). This periodicity property would be lost in an expansion in powers of D . We should 
remark that no restriction has been put on the field configuration, which is completely 
general (may be non Abelian and non stationary) and also no choice of gauge has been 
needed. 
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Perhaps it should be emphasized how exactly the lack of periodicity would break gauge 
invariance. To see this it is sufficient to consider + 1- dimensional fermions in presence of 
an Abelian configuration flu |. The corresponding effective action is of the form (3g(a) with 



a = J (IxqAq. The quantity a is invariant under topologically small gauge transformations, 
Aq = A + doA (A(xq) being a periodic function) but under a large gauge transformation, 
e.g. Aq = Aq + 2nin/j3 (which corresponds to U(xo) = exp(2ninxo/ (3)), a changes by 
an integer multiple of 2ni, so g(a) will not be invariant in general. When g is periodic 
the effective action becomes with Q = exp(— a) and it is invariant under all gauge 

transformations. See Section [IV A| for further remarks. 



D. Relation with the calculation fixing the gauge 



In [ p3| , |rT| the kind of calculation just described was carried out for fermions but fixing 



the gauge through the gauge condition doA = 0. (The idea was that the two operators 
treated not perturbatively, 8q and A , are then commuting.) No loss of generality is actually 
implied by this approach since such a gauge always exists P3"] . However, because it is not 
unique, it is necessary to find all remaining gauge transformations allowed within the A Q - 
stationary gauge, and then check that all of them produce the same result. This was shown 
to be equivalent to the periodicity condition that follows from summing over Matsubara 
frequencies, eq. (p0|). All this is unnecessary in the present approach since the gauge has 
not been fixed. Using gauge invariance, the results obtained within the ^-stationary gauge 
can directly be taken over as follows. When OoAq = 0, the field f2 becomes e _/L4 °, so it is 
only necessary to replace e _/3j4 ° of the calculation in Aq- stationary gauge by Q(x) to obtain 
the result expressed in an arbitrary gauge. Further comments are made in |IV A| . 



E. Expansion in space-time derivatives at finite temperature 

As noted, expanding in powers of Dq breaks periodicity and hence gauge invariance, 
however, in principle nothing prevents from expanding in powers of T>q in eq. (^), namely, 

oo „ 

TT = J2 / d d+1 xtr[ip n (n)XV%Y} , (36) 

n=0 J 

where 

oo 

<p(u,y) = ^Vni^y 71 . (37) 

71=0 

The interest of doing that is, of course, that, at least at lower orders, the result is simpler 
than the full result. This can be regarded as the finite temperature generalization of the 
usual derivative expansion at zero temperature. Recall that D already was restricted to 
appear in commutators only, so this is really an expansion in powers of T>^. As usual, 
higher orders are increasingly ultraviolet convergent. It can be noted that in the Abelian 
and stationary case V vanish identically (on multiplicative operators, such as X and Y in 
eq. (|36"D ) therefore the zeroth order in the above expansion becomes exact. 
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Nevertheless, it should be noted that such an expansion is not as well-defined as for 
instance the expansion in powers of D. The latter is defined from the functional itself, 
since it corresponds to spatial dilatations of the fields. No such transformation is known 
for the expansion in powers of V . So it principle it should be expected that different 
ways of expressing the functional in terms of Q and T>q would yield different expansions, 
only the sum of all orders being unambiguously defined. This can be seen more clearly as 
follows. Recalling that inside (0||0) the operators D 03 and D m are equivalent, the typical 
term considered above, eq. (0), takes the form 

TT = tr(0|<?(A)i,A)2)^nO) (38) 



(with g(x,y) = g(x,y,x)). Using the final form eq. (p3|), it is easily established that it can 
also be written as 

TT = tr((%(A)2, D O i)YX\0), (39) 

because in eq. ( |33"1) all operators are multiplicative and therefore integration by parts and 
the cyclic symmetry can be used. Now, in the frequent case of contributions where X = Y, 
this implies that only the symmetric part of the function g(x,y) is actually contributing. 
However, it is easy to write purely antisymmetric and periodic functions g{x,y) such that 



when used in eq. (p6|) , each order is non vanishing, although of course their full contribution 
vanish when summed to all orders. This particular kind of ambiguity can be fixed by 
imposing a symmetry restriction on g(x,y) before carrying out the expansion in T> . This 



ambiguity is further discussed in Section [IV B 



F. Illustration of the method 



To illustrate the previous manipulations in a practical case, we will consider the quantity 

2" 



C[m, A] 



1, 

4 



: Tr h 



_J \ a F 

-Dj. + m z 2 



(40) 



which will appear later in the study of the scalar field in 2+1 dimensions. In this expres- 
sion F^y = [D^,D V ] and a^ v = §[7^,7^] (where 7 M are Hermitian Dirac matrices in 2+1 
dimensions), m is a c-number. 

C[m, A] is ultraviolet finite and so it is a well-defined and unambiguous quantity. We 
will compute it through second order in an expansion in the number of spatial covariant 
derivatives. Clearly the expansion starts at second order and it is sufficient to retain the 
explicit term in eq. ( ^0[ ) since terms of 0(F^ U ) must contain four or more spatial indices. 

First, the symbols method, eq. ( TfT|) , is applied. Afterwards, taking the trace in Dirac 
space and keeping just terms with two spatial covariant derivatives, produces 



tr 



1 



-(D +p ) 2 + p 2 + m 2 



E 



|0>, 



(41) 



where p 2 = —p 2 , p = 2mn/(3, and E = [D , D] 
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To proceed to the integration over momenta and sum over frequencies we use the trick of 
adding a label 1 or 2 to the operators D Q to indicate their actual position in the expression, 
namely 



-(An + Po) 2 + p 2 + m 2 -(D 02 + Po y + p 2 + m 2 J 1 ' 
:=(0\tr[g(D 01 ,D 02 )E 2 }\0) (42) 

The momentum integration yields 

/ \ _ 1 , / m 2 - (xi +p ) 2 \ . . 

9 ^ X2 >- 8vr (3 ^ (xx+po) 2 - (x 2 +p y g \m 2 ~ (%2+Po) 2 J ' 1 1 

In order to sum over frequencies, it is convenient to reduce the expression to a rational form. 
This is achieved by derivating with respect m and then integrating back (using that g(xi, x 2 ) 
vanishes asm-> oo). Then the identity 

E( J = coth(xi) - coth(x 2 ) , (44) 

n x ' 

can be applied. This produces 

1 1 f°° J J cot Hi( t + x ^)) coth(f(* + ar 2 )) 



g(x 1 ,x 2 ) = -- / dtl—^ + p.p.c. (45) 

lb7T xi - x 2 J m \ 2t + xi-x 2 2t-xi+x 2 J 

Where p.p.c. (which stands for pseudo-parity conjugate) refers to the same expression with 
the replacements Xi — > —x\ and x 2 — > —x 2 . 
Correspondingly, 

C 2 [m,A] = J d 3 xti [tp^V^E 2 ] (46) 

with 

^ y) = W,yj m dt { 2t-ye^-u J - 2t + ye^-u J ) +P - p - C - (4?) 

p.p.c. corresponds to y — > — y and u — > uj^ 1 . This is the final expression which contains all 
contributions to C[m, A] with two spatial Lorentz indices and any number of zeroth indices. 
As expected at finite temperature, it is non-local in time but local in x. Note that C 2 [m, A] 
is an even function of m.f] 



4 This can be shown by noting that <p(u, y; m) — (p(uj, y; — m) is given by the same formula ( (47[) with 
replacement m — > — oo, and then showing that the integrand is convergent and odd as a function 
oft. 
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We can now consider a further expansion in powers of T>q since it respects gauge invari- 
ance. An explicit computation shows that ip{cu,y) is an analytic function of y. At leading 
(zeroth) order in T>q the result from eq. ( |47| ) is 

+ O(X? ). (48) 

Note that this result, unlike the full result in eqs. (f5|[|7]), does not contain an integral over 
the mass (J^° dt). This property holds to all orders in V . 



C 2 [m,A] = t~ ~- / d 3 xir 
ID7T 2m 



I e? m - Q + eP m - Q- 1 ) 



G. The effective action 

In this subsection we summarize some properties of the effective action which will be 
needed later. 

The effective action. The Euclidean effective action is defined as minus the logarithm of 
the partition functional. For non interacting fields it takes the form 

W[M,A] = cTr\og(K(M,D)) = c log Det(fT(M, £>)) , (49) 

where K(M, D) is a differential operator (e.g. the Klein-Gordon operator as in eq. (Q) or 
the Dirac operator) and c some constant. 

Ultraviolet ambiguities. The previous expression needs to be regularized, and a number 
of methods can be used to obtain a renormalized version of it. The key observation is that 
all renormalized versions of the effective action must yield the same ultraviolet finite contri- 
butions and so two such versions can differ at most by a term which is a local polynomial 
in the external fields and their derivatives, with canonical dimension no more than d + 1 
(in d + 1 dimensions). Note that for this discussion a mass m plays the role of an external 
scalar field which happens to take a constant c-number configuration and so, in particular, 
the ambiguity in the renormalized action will depend polynomically on m. This implies that 
any sensible (that is, correctly describing the ultraviolet finite contributions) regularization 
plus renormalization prescription can be used to make the effective action finite; the actual 
effective action describing the physical system at hand will correspond to adding the appro- 
priate local polynomial action to the previous result. Another consequence is that formal 
identities can be applied so long as a violation of them is allowed in the form of a local 
polynomial of dimension d + 1 or less. 

The (-function method. The effective action can be defined through the (^-function 
prescription, namely 

Trlog(if) = -^Tr(K s ) . (50) 
ds s=o 

(An analytical extension in s is understood from sufficiently negative values of s.) When K 
admits a complete set of eigenvectors, Ti{K s ) = J2 n A^, A n being the eigenvalues of K. If 
the calculation has to be made using some expansion it is convenient to use the following 
formula 



r 
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where the path V encloses anti-clockwise the eigenvalues of K but not z = 0. This method is 
practical in actual calculations combined with the symbols method: after applying eq. flllD 
to expand the functional trace, it is straightforward to make an explicit expansion in M and 
D, for instance. This method has been used for fermions in |2T| and for a non-local Dirac 



operator in |22] at zero temperature, and at finite temperature for odd dimensional fermions 



in |23f and even dimensional fermions in 

Anomalies. If the effective action breaks a symmetry of the action there is an anomaly. 
In general the anomaly can be defined as the difference between the effective action of the 
original and the transformed configurations (of the external fields) and by construction is a 
local polynomial. It may happen that the symmetry can be restored by adding an appro- 
priate local polynomial to the effective action. In this case the breaking is an unessential 
anomaly. When not all symmetries can be restored simultaneously the theory presents an es- 
sential anomaly. All symmetries which are implemented through a similarity transformation 
of K(M, D) leave the spectrum invariant, thus, because the (^-function prescription defines 
the determinant of a differential operator by regularizing the product of its eigenvalues, it 
follows that for these symmetries the (- function version of Det(f^) is always free of anoma- 
lies. This applies in particular to vector gauge invariance. (Axial gauge transformations and 
scale transformations, for instance, are not implemented by similarity transformations and 
they can be anomalous.) Therefore the partition functional is always invariant under gauge 
transformations. On the other hand the effective action can change by integer multiples of 
27ri, since the logarithm is a many- valued function. (We are assuming that no zero modes are 
involved. They would induce changes multiples of iir in the effective action.) By continuity, 
this can only happen for topologically large gauge transformations. For a scalar field the 
multivaluation cannot occur since the corresponding ^-function renormalized effective action 
is purely real. For fermions the multivaluation may take place depending on the topological 
numbers of the gauge transformation and the gauge field configuration and this indicates the 
presence of topological pieces in the effective action Such multivaluation is indeed 



found in the explicit calculation for 2 + 1-dimensional fermions of [23] . 

Locality and finite temperature. In previous subsections we have considered operators 
of the form f(M,D) with / one-valued and ultraviolet convergent. Actually, one wants 
to compute the effective action which contains multivaluation and ultraviolet divergences, 
eq. (f£9|). In many expansions (perturbative, derivative, 1/m, etc) higher orders are ultra- 
violet finite and thus they are also free of multivaluation. For those terms all our previous 
considerations hold directly. In particular, we find a remarkable result, namely, that the 
effective action at finite temperature can be written as (cf. eq. (|36])) 

W[M, A] = J2 f d d+1 xtr[ Vn (tyO n ] , (52) 

where tp n are some functions and O n are gauge covariant and local operators constructed 
out of T>^ and M. In this sense the theory at finite temperature is local in the usual sense 
(i.e., the effective action admits an expansion in £> M ) provided that the field Q(x) is regarded 
as local. 

For the ultraviolet divergent terms, some oddities appear which are necessary in order 
to accommodate the existence of anomalies, topological terms and multivaluation, all these 
issues being related. For instance, when the expansion in spatial covariant derivatives is 
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computed for fermions in 2 + 1 dimensions [^3] the functions ip n of lower orders are many- 
valued (a property belonging to the exact result in + 1 dimensions [[l(J and in 2 + 1 
dimension for Abelian and stationary configurations [Q). In addition, negative powers of 
T>q may appear when one goes beyond the Abelian and stationary case. This was handled 
in by introducing the fields A.(x), defined as any solution of the equation V E = VqA.. 
There it was shown that the ambiguity in the definition of A. cancels and all solutions yield 
the same effective action. In terms of the *Alx) and Q(x) the finite temperature effective 
action remains local. 



III. THE SCALAR FIELD IN 2+1 DIMENSIONS 

In what follows we will apply the previous ideas to the computation of the effective action 
of 2 + 1-dimensional scalar particles at finite temperature and density. The Euclidean action 
of the system is that of eq. ([!]). There, the field <f>(x) is a Lorentz scalar and a vector in the 
internal symmetry space which will collectively be referred to as flavor space, with dimension 
Nf. The covariant derivative is defined as = <9 M + A M , where the gauge field A^x) is an 
antihermitian matrix in flavor space. Correspondingly, the gauge transformations U(x) are 
unitary matrices. The mass m is a space-time constant and a real c-number in flavor space. 
The more general case of an arbitrary scalar field M(x) replacing m will not be considered 
here. The effective action is given by eq. @. 

Relevant symmetries of the problem are pseudo-parity and gauge transformations. Pseudo- 
parity corresponds to changing every Lorentz zeroth index, i.e, (xq,x) — > (—Xq,x) and 
A — ► — Aq. Since the spectrum of the Klein-Gordon operator — D 2 + m 2 is unchanged un- 
der this transformation the ^-function regularization prescription provides a pseudo-parity 
preserving effective action. Such Euclidean effective action contains only contributions with 
an even number of Lorentz zeroth indices, it does not contain the Levi-Civita pseudo-tensor, 
and thus it is purely real. Any other renormalization prescription can only produce imag- 
inary contributions which are local polynomials. As noted, the ^-function regularized ef- 
fective action will be strictly gauge invariant since it is real. In fact no essential anomalies 
are present in the case of scalar fields in 2+1 dimensions (scale anomalies are absent in odd 
dimensions p6| ) . 

There is a third symmetry, namely, the transformation m —>■ —m which is trivial for 
scalar particles and again free from anomaly using ^-function regularization. Within other 
renormalization schemes there can appear terms breaking this symmetry but they will be 
removable by adding a local polynomial. In the case of odd-dimensional fermions neither 
pseudo-parity nor the transformation m — > — m are symmetries, however their product 
gives the parity transformation. Parity is a symmetry of the fermionic action but is not a 
similarity transformation of the Dirac operator, so it is not guaranteed to be preserved by 
the (^-function renormalization prescription. As is well-known, parity for odd dimensional 
fermions is in general in conflict with invariance under large gauge transformations and if 
the latter invariance is enforced, parity may present an anomaly, depending on the number 
of flavors f|. 
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A. The 0+1 dimensional model 



The above-mentioned remarks can be illustrated with the 0+1 dimensional version of 



the system. The corresponding effective action has been computed in |L8|. Perhaps the 
simplest way to derive this effective action is by computing the partition function of the 
associated O-dimensional Hamiltonian system in a gauge where Aq is time independent. The 
energy spectrum is then obtained directly from the Klein-Gordon equation written as 



d 4>= (±m + A o )0. 



(53) 



Thus for each flavor there are two single particle levels e± A = ±m + A 0a (a = 1, . . . ,Nf 
labelling the eigenvalues of the matrix A in flavor space). The standard textbook result for 
the partition function of a system of non-interacting bosons then applies 



N 



f oo 



Z s [m,A] = Hm2< 



-(3(o-m+A 0ta )n 



ct=± a=l n=0 

or equivalently 

W s [m, A] = tr log [(l - e -^ m+A ^) (l - e^"" 1 ^)] . 
The trace refers to flavor space. This result can be rewritten as 

W s [m, A] = -/3tr(A ) + T s [m, A] . 



(54) 



(55) 



(56) 



The first term is the 0+1 dimensional Cher n- Simons action which breaks pseudoparity and 
can be removed by a local polynomial counterterm. The second term is (up to a constant) 



T s [m, A] = tr log 



i sinli ( ^(m + Aq] 



sinh 



Ac 



(57) 



This effective action is an even function of m and Aq, so it preserves parity and pseudo-parity. 
It can be written in a manifestly gauge invariant form as 



T s [m,A] =tr log [e^ + e-f 3 ™ 



Q-Q- 1 ] . 



(58) 



As in the case of fermions JT0[, periodicity of the effective action as a function of (3Aq, would 
be lost within a perturbative expansion, i.e. an expansion in powers of Aq. 

Since r s [m, A] enjoys all symmetries of the action it coincides with the (^-function regu- 
larized effective action, up to a constant (since any other local polynomial must be of degree 



one in m or Aq and would break parity). The result in [18] corresponds to r s [m, A] — T s [m, 0] 
in Minkowski space. It is noteworthy that the partition function defined directly from the 
Hamiltonian breaks pseudo-parity (due to the Chern-Simons term in eq. (|56"1)) even if no 
ultraviolet divergences are introduced in 0+1 dimensions in the canonical formalism. 
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B. Computation of the effective action in 2+1 dimensions. Relation to the fermionic 

case 



The effective action of the 2+1 dimensional model cannot be computed in closed form 
for arbitrary space-time and internal symmetry configurations. Our approach will be to 
expand the effective action in the number of spatial covariant derivatives, or equivalently, in 
the number of spatial Lorentz indices. The computation will be carried out through second 
order, that is, we will keep terms with zero or two Lorentz indices. (There are no odd order 
terms in the expansion.) The zeroth Lorentz index dependency is treated exactly since 
this guarantees the preservation of the periodicity condition which is essential for gauge 
invariance. 

The calculation can be done by applying the ^-function regularization prescription with 
the help of the symbols method, as described in |23j , however, it is more economical to use 
the results already established for fermions in that reference. This can be done as follows. 
The effective action for fermions is 



W f [m, A] = -Trflog^D^ + m) 



(59) 



The functional trace is taken in the space of antiperiodic wave functions and includes Dirac 
degrees of freedom, in addition to space-time and flavor degrees of freedom. The gamma 
matrices are hermitian and satisfy 7 M 7^ = 5^ + a^. Actually, there are two inequivalent 
irreducible representations of the Dirac algebra which are distinguished by the label rj = 
±1 in the relation 7^7^ 7 P = ir]e^ up . So if 7^ is one of the representations, — 7^ provides 
another inequivalent representation of the Dirac algebra. The label r\ is attached to the Levi- 
Civita pseudo-tensor and thus a change in rj is equivalent to a pseudo-parity transformation. 
Therefore the fermionic effective action can be split into two components 



W f [m, A] = Wf [m, A] + r\Wj [m, A] 



(60) 



where Wf~ is real and even under pseudo-parity and Wj is imaginary and pseudo-parity 
odd. (Of course, this relation can be violated by adding a local polynomial.) Next, note that 
the formal identity Tr log(AB) = Tr log(A) + Tr \og(B) holds for the functional trace up 
to ultraviolet divergent contributions and so it holds modulo local polynomial terms. This 
implies that 



Wf[m, A] = --Tr f log [(7^ + m)(-7^ + m)\ 



--Tr/log 



-K + m 



{F IH ,= [D ll ,D v )) 



(61) 



So Ws[m, A] (cf. eq. (§)) is closely related to Wf [m, A]. The differences between both 
expressions are i) the Dirac degree of freedom which is absent in the scalar case, ii) the 
different (periodic versus antiperiodic) boundary conditions and iii) the extra term —^cr^F^ 
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which is not present in the Klein-Gordon operator. In addition, a local polynomial action 
can be further added at the end.f] 

The first correction amounts to dividing the fermionic result by the trace of unity in 
Dirac space which is two in the 2+1 dimensional calculation in [p3|l . The second correction 
can also be tackled straightforwardly. The Matsubara frequency (related to do in D ) takes 
discrete values 2iti(n + \)/ f3 for fermions and 2itin/ (3 for bosons, thus the functional trace 
computed in the fermionic Hilbert space is related to the bosonic one after the replacement 
A — ► A — This is equivalent to Q(x) — > —Q(x). (Q(x) has been defined in eq. (f28|).) 

ITT 

W a [m, A] = -W+[m } A — —,A]— C[m } A] . (62) 

P 

The term C[m, A] takes into account the spurious contributions coming from — hcr^F^, 
which have to be removed from the fermionic result. 

This formula can be illustrated in the 0+1 dimensional model, where it reads W s [m, A] = 
— 2WjT[m, Aq — ?j] (note that Dirac space is 1-dimensional in 0+1 dimensions so the factor 
2 is not canceled in this case, and also C[m, A] = 0). The simplest way to obtain the 
0+1 dimensional fermionic effective action is again using the Hamiltonian formalism (fixing 
d Q A = 0). Since there is a single-particle level with energy e a = r\m + A 0t(l for each flavor 
(where rj = ±1 is the Dirac matrix 70 in 0+1 dimensions) it follows that 

W f [m, A] = - log Yl Yl e~^ m+M ' a)n = -tr log [l + e -/»(v*-Mo)] . ( 63 ) 

a=l n=0,l 

The result in [|Hj corresponds (up to a local polynomial) to W/[m, A] — Wf[m, 0] with 77 = 1. 
This version of the effective action does not directly satisfy eq. (|60D, i.e. the pseudo-parity 
transformation Aq —Aq is not equivalent to the transformation rj — > —rj in the previous 
formula. However, subtracting an appropriate //-dependent polynomial^ yields 

W' f [m, A] = -tr log [l + e -^ m+r,Ao) ] = -tr log [l + e^fi"] , (64) 

which does satisfy eq. (|60|). This is the ^-function result p3| . It is readily verified that 
r s [m, A] in eq. fl58[) coincides with — 2Wj + [m, A — ^] plus a polynomial, J dx^tiim). 

In 2+1 dimensions the subtraction C[m, A] can be computed in an expansion in powers 
of \(T^ U F^ V . The term of first order does not contribute (since the trace of a^ v in Dirac 
space vanishes) thus the leading term is that with two powers of F^, namely, the expression 
in eq. fllPP . Because C[m, A] is ultraviolet finite it is free from anomalies, i.e., it enjoys 



Alternatively, one can choose to change m — ► — m instead of 7^ — > — 7^ in the second factor 
in the logarithm in eq. (|6l|), and then relate W s [m, A] to Wf[m, A] + Wf[— m, A]. Up to a local 
polynomial, this procedure is equivalent to the one used in the text. 

6 To wit, 6(—rj) J dxotv(r]m + Aq), which is temperature independent. This is consistent with the 
fact that the finite temperature does not introduce new ultraviolet divergences. 
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all symmetries of the bosonic action. This term has been computed, up to two spatial 
derivatives, in subsection [II F . 

Let us quote the results for the pseudo-parity even component of the effective action 
of fermions in 2+1 dimensions ||23|| . At zeroth order in the number of spatial covariant 
derivatives the result is 



Wf,o[m,A} = -tr(0\ 











-P. 



|0)+p.p.c. (65) 



At second order 



W+ 2 [m,A) 



in 



: tr(0| 



1/2 
2 \p 

,0, 



MtK" 2 ~ An)) - -0o 



2" Vi// (3 T ^2 



I TYl 

m-D 01 )) ( -P02 - An) + y 



G?i tanhf— (t — Dot)) , 

V 0lJJ 2t + D 02 -D 01 J(D 02 -D 01 f 



E 2 



10) 



+X 12 + p.p.c. 



(66) 



In these formulas, p.p.c. means pseudo-parity conjugate, D — > —D , and X\ 2 means the 
same expression exchanging the labels 1 and 2. 
The functions 4> n (z) are given by 



(w) = P n +i(u) 



dz± r-L- (tanh(z) - 1) , n = 0,1,2, 

n! 



Re(w) > 0, (67) 



where the integration path runs parallel to the real positive axis towards +oo. The P n (uj) 
are polynomials of degree n 



Pi(w) =w, 

P 2 (uj) = -UJ 1 — 

K ' 2 6 V 2 



PM 



1 / m 



6 



-uj 



6 V 2 



(68) 



These formulas for 4> n {uj) refer to Re(uj) > 0. When Re (a;) < the property (p n (—uj) 
(— l) n (f) n (uj) can be used. 



C. The effective action 

We can now use eq. fl62"|) to obtain the effective action of the 2 + 1-dimensional scalar field. 
Up to two Lorentz indices, C[m, A] equals C 2 [m, A] which is given in eq. (flop . W^[m, A] 
is given in eqs. (|65| ) and (|66|) . In these latter expressions the functions n (z) are made 
explicit using eq. (]67D , and the variables f2 and P are used instead of _D i and _D 2 (cf. 
subsection |11 (J|) . In addition, we introduce a chemical potential by means of the replacement 
A —>■ A — fi || where \x is a real constant c- number (recall that A is antihermitian). This 
shift is gauge invariant and it is equivalent to Q — > e^fi. 
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The effective action up to two spatial covariant derivatives at finite temperature and 
density is thus 



(69) 
(70) 



W,, [m,A] = J d 3 xtr [Me^fy] , 
W s , 2 [m,A] = [ d 3 xtv[ V2 (e^Q 1: V 02 )E 2 ] 



The functions (po and (p 2 are given by 

j / 2 p+oo 



4vr \3 



Iml 3 + 



dt(t 2 - m 2 )- 



m 



CO 



p.p.c. 



(71) 



(p 2 (v,y) = -tt\ 



1 fl 



-\m\y H — \ —y 



-r / tft 

|m| 

+p.p.c. 



2oo 



m 



e& - uj V 2t - y 



1 1 log 



2 


m 


+ 1/ 


2 


m 


- y 



2uj 



eP(t+y) _ u V 2t + y 



f 2 



(72) 



In these formulas p.p.c. corresponds to y —>■ —y and a; — > a; -1 . 

The formulas (|69|) and (|70|) , expanded in eqs. (|7l|) and (|72|) , constitute the main result of 
this section. They are necessarily complicated looking since an infinite number of Feynman 
graphs (with any number of temporal gauge fields, cf. Section [IV B|) are being added, and 
the effective action is non-local in time at finite temperature. A much simpler formula is 
presented below if only the lowest order is retained in a further expansion in the number of 
temporal covariant derivatives. 

Some remarks are in order. In writing the formula we have already used the fact that 
the effective action is an even function of m (since this is already true for Wt[m, A]). The 
condition Re (a?) > in eq. (|67|), implies that the previous formulas refer to the physically 
relevant case < \m\ only. All polynomial contributions introduced by P n (oj) in eq. (|6~7D 
combine in such a way that all the polynomial dependence on Dqi cancels and the result is 
a periodic function of Dq\. This cancellation requires the explicit and as well as the p.p.c. 
terms to hold and it is a non-trivial check of the formulas. (Of course, the periodicity holds 
also for the fermionic effective actions in the pseudo-parity even sector, W^[m, A].) The 
argument of the logarithm in (p 2 (uo,y) is to be taken in the interval (— 7r,7r) and t runs on 
the real positive axis. (Note that y represents V and thus it is purely imaginary.) 



D. Expansion in space-time derivatives at finite temperature 

A simpler expression for W Sj2 [m, A] is obtained retaining only the leading order in an 
expansion in T> . This corresponds to expand in powers of y in the function (p 2 (u,y). As 
noted this expansion does not break any symmetry and is a natural one in the present 
context. This produces 

W a , 2 [m,A] = [ d 3 xtr [tp 2fl {e^Q)E 2 ] +O(V ), (73) 
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with 



1 1 e 2/3m 
967T m 

1 1 d 
967r m dm 



- LO 
(ePm _ u 



2/3me 



(3m 



LO 



+ p.p.c. 



LO 



m- 



,/3m 



LO 



+ p.p.c. 



(74) 



Note that <£> 2 ,o( w ) is an even function of m. This formula is manifestly invariant under all 
symmetries of the action. Two non-trivial checks of the calculation are i) that negative 
powers of y have canceled and ii) that y^.oO^) no longer contains parametric integrals (i.e., 
J, | dt). As noted, this expression is exact for 14^2 A] in the Abelian and stationary case. 



E. Zero temperature limit 



The zero temperature limit of eqs. (|7lf) and (ff2 ) is straightforward: because < \m\ 
and Q is unitary, the term always dominates, thus the limit j3 — > oo just removes the 
terms with Jj^, dt and the dependence in \i is also canceled (we are assuming \i fixed as 
T —> 0). That is 
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(75) 
(76) 



At zero temperature Lorentz invariance is a symmetry of the action. Such symmetry is 
not obvious in W Sj 2[m, A] since it sums all orders in T>q but not in T>i, i = 1,2. However, 
considering an expansion in inverse powers of m (a Lorentz invariant expansion) gives 

W s>2 [m, A] = ~^J^J d3xtT [E 2 ] + 0(m~ 3 ) . (77) 
3D. This expression admits a unique Lorentz invariant 



The same result follows from eq. 
completion, namely 



W s [m, A] = - d 3 xtr 



— m + T-iFt, 

6n 96n \m\ 



+ 0(m~ 



(78) 



This formula comes from a direct application of the results in [g7]] thereby being a check for 
our formulas. At next order in 1/m several Lorentz invariant operators of dimension 6 can 
appear and eq. (f76| ) puts a constraint on their coefficients. Note that the last two formulas 
hold also at finite temperature since the temperature-dependent corrections are 0(e _/3 ' m ') 
(see also Section [IV I} ) 



F. The partition function 

The effective action at finite temperature and density is directly related to the grand- 
canonical potential, namely W s [m,A; (3, (/] = f3VL(f3,n). (Note that \i introduced by the 
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replacement do — » <%— /i couples to the charge and not to the number of particles which is not 
conserved in the relativistic case [p8]|). However, strictly speaking, a system at equilibrium 
with temperature T and chemical potential /i can only be stationary. In addition, the 
physical effect of a (negative) constant external scalar potential is indistinguishable from 
that of a (positive) chemical potential, since both add to the energy for positive charges and 
subtract for negative ones, thus for the partition function Aq should not be Wick rotated, 
and Aq is real instead of imaginary. All expressions depend only on the combination Aq — fi. 
(It would be algebraically inconsistent not to rotate A to its Euclidean version in the general 
case, but not within the subset of stationary configurations.) Note that the fact that A is 
real or imaginary does not affect the functional form; in any case the functional derivative 
of the effective action with respect to Aq yields the charge density and the derivative with 
respect to /i yields the total charge: 

In view of the relation between partition function and effective action, it follows that 
W S fl describes a two-dimensional relativistic ideal gas in presence of an almost space-time 
constant scalar potential Aq. An explicit calculation of the charge density using W s q in 
eq. (|TID and integrating by parts yields 

jy poo (3 8 \ 

P» = jij l d«^coth(|(t-/i))-coth(|(t + /i))J , (80) 

Nf = tr(l) is the number of flavors. We are assuming that fi couples equally to all flavors 
and we have dropped Aq since it can be recovered from /i. This formula can rewritten in 
the standard form |28] 



i d 2 k f N f Nf * 



where u(k) = Vk 2 + m 2 . 

Likewise, the term W s ,2 adds a contribution to the density coupled to Aq. Let us con- 
sider the Abelian case (in addition to stationary), then T>q vanishes and eq. ( [73]) becomes 
exact. There will be two contributions to the charge density, one coming from the explicit 
dependence on Aq and another though the dependence in E. The latter is a total derivative 
and does not contribute to the total charge. Both contribution can be combined to give 



j x ) = _^LlA( m 
967r m dm 



- ^cosech 2 (^(m - fx + A ))E 2 



+2 coth(^(m -n + A )) V.E 



p.p.c. (82) 



(Note that, A , \i and E are pseudoparity odd.) At zero temperature (assuming \m\ > 
\n — Ao|, or else at finite temperature but large mass) this simplifies to 
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which also follows from eq. (|77Q . 

It is also interesting to note the relation of our results with other formulas using polylog- 
arithms.[] The functions 4> n (u) introduced in [^3] and subsection [111 U| are directly related 
to polylogarithms, namely 



,{u) = P n+1 (oo) - (-2)-"Li n+1 (-e- 2w ) 



(84) 

This relation is easily established from eq. (|67| ) by noting that it verifies the following defining 
properties of Li n (z): 



hin(z) 



z^-lX n {z) = Li n _i(z) , Li„(0) = . 



In this notation, W s ,q in eqs. flBTf ) and (fn]) becomes 



W a ,o[ra,A] = ~7~y d 3 xtr 
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2T 3 {P\m\Li 2 {z)+U 3 {z)) 



+ p.p.c. 



(85) 
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where z = e ^' m ' is the fugacity, and the corresponding density becomes 



N f T 2 
2vr 



|m|Lii(z) + Li 2 (z)) - p.p.c. 



(87) 



in agreement with [30 



As a final comment related to the partition function, we note that the relation eq. ( p2|) 
can also be written using the chemical potential instead of the scalar potential Aq. That is, if 
£lb{P, A 4 ) and Qf(/3, //) represent the grand-canonical potentials of a system of non interacting 
particles in presence of external fields, treated as bosons or fermions respectively, then 



Z7T 



(88) 



The minus sign comes because the functional integral with Grassman variables gives the 
determinant of the quadratic form instead of the inverse determinant. The shift fi — * fi+ ^ 

accounts for the different boundary conditions. Of course a shift \x — > \i + must leave 
the partition function invariant, since \x is coupled to an integer-quantized charge. This is 
another manifestation of gauge invariance. For interacting particles eq. (|88|) can be extended 
using the well-known prescription of adding a minus sign for each particle loop (eq. ( |88"|) 
corresponds to the particular case of one-loop). 



7 The polylogarithms are defined as [29| 



oo u 



Li n (z) = ^ — 



fc=l 



8 An equivalent procedure would be to compute the grand-canonical potential for n replicas of the 
particles and then set n to —1 at the end. That is, in the notation of [MJ the Hamiltonian becomes 
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IV. FURTHER GENERAL REMARKS ON THE METHOD 



A. Remarks on gauge invariance 

As noted in the introduction, topologically large gauge transformations have played a 
prominent role in the development of this field by putting severe constraints on the allowable 
forms of the effective action functional. On the other hand it seems that in our present 
approach such a role is played instead by the periodicity constraint which prompts the 
appearance of the Polyakov loop Q(x). In this subsection we will make some remarks to try 
to clarify the relation between both concepts. 

Let G be the gauge group, and let T be the relevant homotopy group controlling the exis- 
tence of topologically large gauge transformations at finite temperature. In + 1-dimensions 
T = 7Ti((j), therefore there are large gauge transformations for the Abelian group U(l) but 
not for simply connected groups such as SU(n). In d + 1-dimensions (with d different from 
zero), let us assume for this discussion that the spatial boundary conditions are such that 
the space becomes effectively compactified to a sphere. Hence V contains classes of map- 
pings from the space-time manifold S 1 x S d into the gauge group G, where the factor S 1 
corresponds to the compactified Euclidean time and S d to the <i-dimensional space. (Note 
that, unless otherwise stated, we will require all functions to be continuous on the space-time 



manifold, and in particular, periodic as a function of time.) Theorem 4.4 of ]52] then implies 
thatf] 

r/7r d+1 (G) = 7r 1 (G)x7r d (G). (90) 

Let us consider the case d = 2. Because the group ^(G) is always trivial for any (compact) 
Lie group G, this simplifies to r/7T3(G) = vri(G). Two typical cases are 

i) G = U(l). In this case 7r 3 is trivial and iri is Z, so T = Z. There are non-trivial gauge 
transformations which wind n times around the S 1 factor of the space-time. They can 
be realized by space-independent, but time-dependent, gauge transformations. 

ii) G = SU(n) (n > 2). In this case 7i"i is trivial but 7r 3 = Z and so once again T = Z. In 
this case the corresponding large gauge transformations must be space-time dependent. 

Note that time-independent gauge transformations are controlled by the homotopy group 
7r 2 (G) which is trivial, and so they are always topologically small in two spatial dimensions. 

Let us now turn to the point of view used in this work. The periodicity constraint 
refers to the fact that a gauge invariant functional must depend on Q(x) and not just on 
\og(Q(x)). The cleanest way to formalize this is by working on the gauge doA = 0, in 



H = Y h <*p Y a o 



T a/3a+ Y ' 
a,/3,7,<5 
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(89) 



0"l,cr 2 = 



9 When d = the theorem is consistent with V = tti(G) x it\{G) which follows from S° = {1, — 1}. 
Physically, we want the spatial manifold at d = to be just {1} and so Y = tt\{G). 
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which Q = exp(— (3Aq(x)). Thus both Q and Aq are time-independent in this gauge. Let 
us remark that taking A to be time-independent is not merely a restriction on the set of 
possible gauge field configurations, it is a choice of gauge in the sense that every gauge 



field configuration admits a gauge transformed configuration which is A -stationary ||23|| . In 
this gauge the periodicity constraint expresses that a gauge invariant functional must be 
a periodic functional of Aq. To be concrete consider a generic expression of the form (cf. 
eq. (0)) 

F[M,A] = J2 I d d+l xtr[g n (A )O n } , (91) 

where O n are gauge covariant local operators, then gauge invariance requires the functions 
g n (z) to be periodic with period 2m/ (3. The necessity of this requirement follows immediately 
from considering the following class of gauge transformations 

U(x) = exp(x A(a;)) , (92) 

where the time-independent function A(x) takes values on the Lie algebra of G and is 
restricted by the following conditions 

L4 (aO, A(a;)] = , exp(/3A(aj)) = 1 . (93) 

The second condition means that the eigenvalues of A(cc) are of the form Xj = 2mrij/[3, 
for integer rij (such integers are ^-independent by continuity) and it ensures that the corre- 
sponding U (x) is periodic in the temporal direction. Under such a gauge transformation 

A T {x) = A {x)+A(x), (94) 

i.e. the spectrum of A is shifted by multiples of 2m/ f3 and the functions g n (z) must be 
periodic. 

Before proceeding, an important point should be noted regarding the approach used 
in this work. Namely, at the price of working with asymptotic expansions, we can afford 
to derive formulas which are "universal" in the sense that no restriction is put on the 
algebraic properties of the internal space, in particular, the formulas must hold for any 
gauge group. This means for instance that the functions g n (z) above are the same for all 
theories and configurations. The requirement of universality puts stronger constraints on 
the functionals that cannot be appreciated when working with concrete theories only (for 
instance, in particular theories some of the operators O n can vanish identically and so the 
corresponding function g n (z) does not play a role). 

The transformations introduced in eq. (|9~2"D, subjected to the conditions eq. (|9~3"D, have 
been named discrete transformations associated to A (x) in |23 since in general they form a 
discrete set due to the condition Xj = 2mnj/(3. Note that this set depends on the particular 
time-independent field A , through the condition [Aq, A] = 0. Clearly, these transformations 
leave invariant the gauge condition doA = 0. Likewise, the gauge condition is also preserved 



by time-independent gauge transformations. In p3[ it is proven that these two kinds of 
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transformations are the most general ones which preserve the Ao-stationary condition. This 
means that within this gauge a functional such as T[M, A] above is gauge invariant if and only 
if the functions g n (z) are periodic. Then it can be written in a manifestly gauge invariant 
form (without gauge fixing) as in the right-hand side of eq. (|5^) with ip n (e~ l3z ) = g n (z). 
Note that invariance under time-independent gauge transformations do not impose further 
constraints on the g n (z). 

The previous discussion suggests a comparison between large gauge transformations, on 
the one hand, and discrete transformations, on the other. Two questions pose themselves at 
this point. Are the discrete transformations large in the topological sense? Is it necessary 
to rely on large gauge invariance in order to arrive to an Ao-stationary gauge? 

The first question can be answered as follows p3|| . For a multiply-connected group such 
as U(l), the non-trivial discrete transformations are topologically large since they loop once 
or more on the temporal circle S 1 . On the other hand, for a simply-connected group such as 
SU(n) (in more than one space-time dimension) the discrete transformations associated to 
some gauge configuration may be large or small, depending on A (x). For instance if A (x) 
is everywhere diagonal, A(cc) is also diagonal and in fact a constant. In this case the discrete 
transformation describes a single loop on the gauge group for all x and it is homotopically 
trivial. In general, however, discrete transformations can be topologically large. It might 
seem that the form of the discrete gauge transformations in eq. (|92|) factorizes time and 
space and so it is always classified by the homotopy group H\{G) being always small for a 
simply connected group. This is true when A(x) is constant or homotopic to a constant, but 
in general this is not the case. The reason is that although A(x) is a map from S 2 into the 
Lie algebra of G (a vector space, and thus contractile), it cannot be contracted to a point 
since the spectrum of A(x) is constrained to be in ^pZ. An explicit SU(2) example in 2 + 1 
dimensions is provided in [E3f, namely, 



f 27mxo \ 

U(x) = exp I — — — irx I , n eZ, (95) 

where r are the Pauli matrices and x lies on the unit sphere S 2 in M 3 . U(x) covers SU(2) 
2n times and thus it is homotopically large for non-vanishing n. 

Regarding the second question, whether a given gauge configuration can be brought to a 
Ao-stationary gauge using only small transformations, it also depends on the group and the 
configuration |p3 |. For the group U(l), any gauge configuration is in the same homotopy 



class as one which is A -stationary. For a simply-connected group such as SU(n), it depends 
on the initial A (a;).[^] 



10 This is the generic case which we will assume. It holds whenever exp(/3Ao(a;)) is either nowhere 
degenerated or at least the regions of degeneracy are sufficiently small that a unique eigenbasis 
can be selected (up to normalization) by continuity |23|| . In this case Aq(x) is also nowhere 
degenerated and thus A (a;) is completely determined by its eigenvalues. So generically the discrete 
transformations form a discrete set. 

11 This can be seen as follows. Let Aq(x) be some ^o-stationary configuration for which all its 
discrete transformations are small, and let us further assume that time-independent transformations 
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An apparent paradox arises here. As emphasized in |IT|,13| , although perturbation the- 
ory breaks large gauge invariance, it respects invariance under small transformations, or 
equivalent, under infinitesimal ones. (Actually, this is strictly correct for Abelian theo- 
ries only. In non- Abelian theories infinitesimal gauge transformations mix different orders, 
however, the mixing is mild since it only involves finite sets of Feynman graphs.) On the 
other hand, we have just seen that expanding a functional such as T[M, A] in powers of 
Aq (which is a perturbative expansion) destroys periodicity and thus gauge invariance un- 
der discrete transformations. This looks paradoxical since the breaking occurs even if the 
discrete transformations are topologically small, and no breaking was expected in this case. 

The resolution comes from the observation that the right-hand side of eq. ([)!]) refers solely 
to configurations in the gauge 8qAq = 0. Within this gauge the only allowed infinitesimal 
transformations are the time-independent ones, for which no breaking occurs. Non trivial 
discrete gauge transformations, even small ones, cannot be reached continuously within 
the A -stationary gauge. In this sense they are always topologically large. It should be 
realized that the concept of homotopically trivial is a relative one. A transformation which is 
topologically large within the gauge group G can become small if G is regarded as a subgroup 
of a larger group G' and deformations within G' are allowed. Conversely, a small discrete 
gauge transformation becomes homotopically non-trivial if one insists on preserving the 
gauge condition 8qAq = 0. The usefulness of the concept of small and large transformations 
remains: perturbation theory provides a functional valid in the region of small fields, this 
region is preserved by infinitesimal transformations and thus the perturbative functional 
must be small gauge invariant. Large transformations, on the other hand, necessarily move 
the gauge configuration away from the perturbative region and thus the response of the 
perturbative functional under large transformations is not trustworthy. As a consequence 
invariance under large transformations provides useful non-perturbative information and 
puts non-trivial constraints on the functional. This holds whether the transformations are 
large from the point of view of G or from the point of view of the submanifold of Ao-stationary 
configurations. 

For another argument, we can recall our previous remark that a functional such as 
r[M, A] in eq. ( pT|) must hold for all theories at the same time. It is not surprising to 
find a breaking of gauge invariance under discrete transformations in a perturbative expan- 
sion, when such transformations happen to be large, and conclude that non-trivial non- 
perturbative conditions, namely, periodicity, are required to avoid the breaking. However, 
as we have emphasized, the property of being topologically large or small depends on the 
group, whereas the formula must hold in all cases, and so periodicity must follow in all cases 
too. 

In this subsection we have considered a gauge fixing condition in order to deal with 
the quantity logf2 in a simple way. We must recall, however, that the expressions derived 
with the method studied in this work are all fully gauge invariant, provided some regularity 



are also small (for instance d = 2). It follows that all ^-stationary configurations related to the 
previous one by a gauge transformation are in the same homotopy class. Thus, if Aq (x) is a 
gauge transformed configuration with U(x) large, no small transformation will bring it to the 
Ao-stationary gauge. 
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conditions are met, since they depend on f2. This does not mean that large or discrete gauge 
transformations play no role whatsoever. This is because the regularity conditions fail for 
the ultraviolet divergent pieces of the effective action. This translates into the fact that the 
functions f n (fl) (cf. eq. fl52j) ), can be many- valued (although Q = needs not be one of the 
branching points). A typical example is the effective action of + 1 dimensional fermions, 
eq. fl64|). Further let us consider the Abelian U(l) case, which admits large transformations 
and they coincide with the non-trivial discrete transformations in the Ao-stationary gauge. 
Choosing r] = +1 the branching point is at Q = — exp(/3m). Under a large transformation, 
the argument may go into a different Riemann sheet (depending on the sign of m), yet the 
functional is such that it changes, at most, by an integer multiple of 2m, reflecting the fact 
that the partition functional is one- valued. This is a general property which puts constraints 
on the functions ip n (Q). 

Perhaps this is a good place to remark that the gauge invariance implied by the use of 
the ^-function prescription, eq. (|50D, has two levels f2"3|| . One the one hand, there is the 
somewhat trivial gauge invariance implied by the fact that the regularization depends solely 
on the spectrum. Since the spectrum of the operator K is left unchanged by gauge (or more 
generally similarity) transformations, gauge invariance follows. However, the definition of the 
^-function introduces a branch cut in the manifold of operators K, each singular operators 
being a branching point on such a manifold. On a given Riemann sheet, the effective action 
functional may display a jump discontinuity along the branch cut. Because the determinant 
of K, being the regularized product of eigenvalues, is a smooth functional, it follows that 
the jump must be an integer multiple of 2m. This is a tighter constraint on top of the 
trivial gauge invariance noted above. For instance, we have noted in the introduction that 
perturbation theory for fermions at finite temperature yields a Chern-Simons term which 
is renormalized by a temperature dependent coefficient. Under large gauge transformations 
this would introduce an unacceptable change in the effective action by a quantity which is 
not a multiple of 2ni. It would be tempting to "solve" the problem by simply replacing the 
Chern-Simons term by a suitable gauge invariant version of it, namely, the //-invariant. This 
prescription restores gauge invariance but introduces jumps which again are proportional to 
the temperature dependent coefficient and thus it can be ruled out. The exact result known 
in particular but non-trivial cases [|IOH T2ll shows that this is not the correct mechanism and 
that the determinant is a continuous functional, with no jumps. 



B. Feynman graphs and large gauge invariance 



We have already noted in subsection [111 C] that expressions such as those in eqs. ( |5§| , |T0D 
involve an infinite number of Feynman graphs. It seems interesting to understand which 
Feynman graphs are being added and gain some insight on how preservation of full gauge 
invariance is related to this. 



To this end, we will first consider the simpler case of + 1-dimensional fermions |L0 
The corresponding exact effective action is given in eq. ( |6^) and that formula holds for 
arbitrary gauge fields which need not be Abelian nor stationary. As noted, when the gauge 
group is not simply connected it supports large gauge transformations, which augment the 
value of the effective action by 2mkn, n,k G Z, where n is the winding number of the gauge 
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transformation and k depends on the gauge group and the sign of the mass (in higher spatial 
dimensions it may also depend on the homotopic class of the gauge field configuration). 

Because the + 1-dimensional formula is exact, all Feynman graphs are included in 
this case. This suggests that all graphs are required in order to reconstruct the Polyakov 
loop Q appearing in the formula. We conclude that large gauge invariance does not act 
selecting a certain subset of graphs. The same conclusion is expected to hold in higher 
dimensional formulas since the dependence on Q there is qualitatively similar to that of the 
one-dimensional case. 

To further discuss this point let us restrict ourselves to the case of an Abelian gauge 
group. In this case large gauge transformations correspond to discrete shifts of a = — log Q = 
J (IxqAq and 

Wf (a + 2vrm) = W f (a) + 2mkn , (96) 

where the integer constant k is known. This equation contains all the information on large 
gauge transformations, and it is completely equivalent to the statement 

W f {a) = P{a) + ka, P(a + 2mn) = P(a) . (97) 

Therefore, large gauge invariance is equivalent to the strict periodicity of the function P(a) = 
Wf(a) — ka. Feynman graphs correspond to expand in powers of a, and looking for large 
gauge invariance in terms of Feynman graphs corresponds to detect periodicity of a function 
from its Taylor expansion. This seems to be a difficult task. 

A related issue is studying to what extent large gauge invariance determines the effective 
action. In the previous + 1-dimensional Abelian case we have seen that to comply with 
gauge invariance, P(a) must be periodic, i.e., 

P(a) = ^c n e na , (98) 

for some (infinite number of) coefficients c n . No further information can be extracted from 
gauge invariance, and the coefficients c n are not determined. In order to achieve further 
restrictions on the function P(a), more information has to be provided. If, for instance, 
one knows that the partition function Z(a) = exp(— Wf(a)) (a periodic function) contains 
only a finite number of periodic modes, the corresponding Fourier coefficients can then be 
determined from a few perturbative terms. This is actually the case in + 1-dimensions 
T9fl , and can be traced back to the fact that the corresponding Hamiltonian contains a finite 



number of states, namely, the vacuum or the 1-fermion state (cf. Section 11 1B|) . In higher 
dimensions, besides the number of fermions, there is a momentum quantum number and 
a corresponding kinetic energy contributing to the eigenvalues of the Hamiltonian, thus in 
general the partition function will contain all kinds of Fourier modes. To see how this works, 
it is sufficient to consider fermions in the Abelian case with Ai = and Aq a space-time 
constant (such A cannot be gauged away at finite temperature). Let e° k denote the single- 
particle levels of the Hamiltonian when A is set to zero (that is, the kinetic energy only) 
and let be the levels when A is switched on. Clearly, e k = e° k + A , and the label k is 
related to the momentum of the fermion. The partition function is 
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Z[m, a] = + e~P ek ) = + fle"^) , U = e~ 
k k 



(99) 



Z[m,a\ will be periodic in a but with many periodic modes, unless d = 0. When d = 2, 
integration over the label k yields a particular case of the formula for W/j !0 [ m > A] i n ec i- (|65D- 

In more than one dimension we have proposed an expansion in the number of spatial 
covariant derivatives. The zeroth order (e.g. eq. fl69])), which contains no spatial indices, 
corresponds to the sum of all Feynman graphs with the spatial momenta pi (of the external 
gauge fields) set to zero and no legs with spatial Ai, but any number of A Q external legs and 
the full dependence in the frequency Pq. There are no odd-order terms in the expansion. 
The second order (e.g. eq. ([7UD) corresponds to all graphs with two Ai or pi, that is, i) 
graphs with Pi = and two spatial gauge fields, plus ii) graphs with one Ai and pi kept up 
to first order in a Taylor expansion of the Green functions, plus iii) graphs with no legs Ai 
and pi kept up to second order in a Taylor expansion around zero spatial momentum. Note 
that the method assumes that the space has topology M d . The contribution of higher orders 
follows a similar pattern. Eventually, all contributions, all Feynman graphs, are added up. 

The previous conclusions follow from inspection of the formulas or else from making 
the expansion by introducing a bookkeeping parameter in order to count spatial indices, 
as explained in Section [II C| . For the simpler formulas obtained by further expanding in 
powers of V , (e.g. eq. (^)), we have noted that this expansion does not seem to follow 
from insertion of a bookkeeping parameter and so different expansions can be obtained, 
all of them being equivalent when added to all orders. The problem and its interpretation 
in terms of Feynman graphs can be seen in a particularly simple case: let us assume that 
the gauge group is Abelian, that Aq is a space-time constant and Ai are space-independent 
although time-dependent. In this case the Green functions depend only on the frequencies 
K n = 27rin/f3 of the fields Af. 

W = f (Q) + J2 Kn)A,nAi,-n + ' ' " , A(x ) = 4 2 ^ A ^ " ( 10 °) 

n n 

In addition, T>q is equivalent to do and so an expansion in T>q is just a Taylor expansion 
in powers of K n , to be made on top of the expansion in powers of A4. This gives the 
interpretation in terms of Feynman graphs. Because n n is a discrete variable, the function f2 
(and similarly for higher orders) is only well-defined at those discrete values. The ambiguity 
comes when it is smoothly extended to continuous values of K n in order to carry out the 
Taylor expansion. Presumably this ambiguity can be removed by using the choice suggested 
by Carlson's theorem. 



C. Large mass expansions 



Large mass expansions of the effective action can be also considered, as done in [13] for 
fermions using the heat-kernel technique. Inspection of our formulas in [111 C| and [llllt| , show 
that as the mass becomes large the temperature dependence is exponentially suppressed, 
being 0(e _/3 ' m '), thus the large mass expansion is an asymptotic expansion with temperature 
independent coefficients. This is consistent with the fact that the coefficients are local 
operators independent of the global topology of the space-time manifold [13|. Therefore, in 
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order to carry out a large mass expansion one can simplify and start with a zero temperature 
theory The simplification of working at zero temperature is enormous to the point that 
this problem can be considered a solved one. There is a very large body of work on this 
subject, both for bosons and for fermions, largely summarized in |3^| and references therein. 



Large mass expansions for fermions with arbitrary Dirac operators and arbitrary dimension, 
computed along the lines of the method discussed in this paper for finite temperature, can be 
found in |2T[ . The large mass expansion can also be used as a check of the finite temperature 
formulas. In |IIIE| we have already noted that eq. (f78[) is consistent with the result more 
straightforwardly obtained from the zero temperature method presented in |J27[| . Likewise, 
for 2 + 1-dimensional fermions, and starting from the full finite temperature calculation, one 
finds 



W + = [ d 3 xtr(F*) + o( -At ] , (101) 

48vr \m\ J v \m 3 J ' K J 

W~ = V ae(-am)W cs - [ d^MFp^M + o(\ 

Sit 12m 2 J ym* 

for the pseudo-parity even and odd components, respectively. (In these formulas Wcs is the 
Chern-Simons term, and e denote the step and sign functions,respectively, r] = ±1 depends 
on the irreducible representation of the Dirac gamma matrices taken and a = ±1 distinguish 
the two possible (^-function definitions of the effective action, depending on the branch cut in 
the function z s .) These results also derive more directly from the zero temperature formulas 
in [^TJ. The term with F^ u in W + is that with if 4 in eq. (4.11) of | I3]| , already noted there, 



whereas the term with Fp tl T> a Fp u is that with P 5 in eq. (4.12) of the same reference. 



V. SUMMARY AND CONCLUSIONS 

Our findings can be summarized as follows. 

1. The one- loop effective action at finite temperature and density, for bosonic or fermionic 
particles in presence of background fields (both gauge and non-gauge) with arbitrary internal 
symmetry group and arbitrary space-time dependence, can be written as a sum (an asymp- 
totic series in general) of terms ordered by the number of spatial Lorentz indices. Each term 
is well-defined from the effective action functional itself and is separately gauge invariant 
under all gauge transformations. 

2. These terms are amenable to explicit computation using a combination of (^-function 
and symbols method. We have shown that this kind of calculations can be carried out 
preserving full gauge invariance throughout, without assuming particular internal symmetry 
groups or special space-time configurations for the background fields. The same arguments 
show that previous calculations done fixing the gauge through the condition doA = can 
be repeated lifting this condition, and this is equivalent to rewrite the final original result 
in a manifestly gauge invariant way. 

3. A further gauge invariant expansion can be taken in the number of temporal covariant 
derivatives in the adjoint representation. Within this expansion, all ultraviolet finite terms 
and more generally, all terms not related to essential anomalies, can be written as a sum of 
gauge invariant local operators (i.e., constructed with = [£> M , ] and M) times a function 
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of the field Q(x), which also transforms locally under gauge transformations. For those 
terms containing anomalies, topological pieces and multivaluation, the effective action still 
looks local in terms of a suitable gauge covariant version of A(x), in addition to M, and 

a 

4. The method is explicitly applied to the problem of relativistic scalar particles in 2 + 1 
dimensions. The corresponding effective action is computed up to terms with two spatial 
Lorentz indices. The result is checked against the known result at zero temperature and also 
the known partition function of a relativistic Bose gas. The corrections to the density are 
also computed. Finally a simple rule is noted relating the bosonic and fermionic versions of 
the grand-canonical potentials of ideal or interacting systems. 
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